Specific heat and orbital moment of CoO from first-principles atomistic calculations 
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We have for the first time calculated low-energy electronic structure both in paramagnetic and 
magnetic state as well as zero-temperature properties and thermodynamics. We consistently de- 
scribed magnetic properties of CoO in agreement with its insulating ground state. The orbital 
moment of 1.42 /j,b gives 35% contribution to the total moment of 4.04 fis at T =0 K. We have 
calculated from this low-energy electronic structure the temperature dependence of the specific heat 
being in nice agreement with experimental data. In our approach CoO is an insulator independently 
on distortions and the magnetic order. 

PACS numbers: 75.40.Cx; 71.10.-w 
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I. INTRODUCTION 

Understanding CoO and other 3d monooxides is impor- 
tant for the development of the solid-state theory. Their 
magnetism and the insulating ground state is a subject of 
a fundamental controversy by more than 70 years. Most 
models basing on the conventional band theory predict 
FeO, CoO and NiO, having incomplete 3d shell, to be a 
metal, whereas they are observed to be very good insu- 
lators, if stoichiometric. 

Despite of 70 years of intense studies the problem 
of the magnetism and the insulating ground-state of 
open-shell transition metal oxides is still under debate 
0, 0, i, 1, i, H, 0, E i M El, 0, El . We have largely 
solved this problem working in the purely ionic picture 
(we call this approach Quantum Atomistic Solid State 
Theory QUASST). In this paper we present theoretical 
description of the magnetism and low-energy electronic 
structure of CoO and the temperature dependence of the 
specific heat resulting from this structure. We earlier 
have described NiO [ll M, \M_and FeO [13 followed 
by earlier description of FeBr9 [l8l| and LaCoOs [lj| . In 
Refs 14-16 we have described NiO within the localized 
atomistic paradigm with only three parameters: the oc- 
tupolar crystal-field parameter B4 =+21 K (lODq =1.086 
eV), the spin-orbit coupling A s -o = -480 K and a small 
trigonal distortion B% =+50 K [li HE |j"6fl. These pa- 
rameters are related to the lowest 21-fold degenerated 3 F 
term. A fourth parameter, the molecular- field coefficient 
ndd of —200 T I [is describes the magnetic state. This 
ndd value yields TV of 525 K in good agreement with the 
experimental observation for NiO. 

A direct motivation for this paper is presentation by 
Kant et ai. [2(| of new measurements of the specific heat 
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FIG. 1: Calculated splitting of the twelve lowest states of the 
4 T\ g subterm for the Co 2+ ion (3d 7 configuration) both in the 
paramagnetic and magnetic state below Tjv of 291 K. The 
used parameters are specified in the figure. The moment of 
the ground state of 4.04 /is is built up from the spin moment 
of 2.62 fiB and the orbital moment of 1.42 hb\ the splitting 
in the paramagnetic state of the highest states is due to the 
spin-orbit coupling - an effect of the tetragonal distortion is 
very small causing only the thicking of the second and third 
levels. 

of CoO with an attempt to describe properties of CoO 
within an ionic model. Also works of Haverkort et al. 
[2l1 | and Wdowik and Parlinski [22j from the last year 
contribute to studies of CoO. 



II. THEORETICAL OUTLINE AND RESULTS 

We have performed detailed calculations of the low- 
energy electronic structure related to the 28-fold degen- 
erated 4 F term [23[ both in the paramagnetic and mag- 
netic state. We took into account the experimental crys- 
tallographic structure - the NaCl-type structure with 
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FIG. 2: The calculated temperature dependence of the spe- 
cific heat of CoO (thick line) contains the magnetic/Schottky 
contribution, the lattice and dilation contributions; thin line 
- the 3d contribution Cd(T) containing both magnetic and 
Schottky contributions. The dashed line shows the lattice 
contribution with single 9d of 600 K. This lattice contribu- 
tion approaches a value of 6 R, R is a gas constant, for high 
temperatures according to the Dulong-Petit law. The dilation 
contribution is approximated by a linear term with tempera- 
ture with -yd =3.4 mJ/K 2 mol. Points are experimental data: 
o - taken from Ref. j2jj, ♦ - taken from Ref. 



a=424 pm and a tetragonal distortion with c/a=0.988. 
We have calculated the strenght of the octahedral crys- 
tal field 10L>t7=1.19 eV (S 4 = -23 K = -2.0 meV - the 
minus sign is related with the octahedral position of the 
Co 2+ ion in the O 2- octahedron; the calculated octupo- 
lar charge moment of surrounding charges at the Co site 
Aj = +267 Ka^> 4 , a^ is the Bohr radius; (r d ) for the 
Co 2 + ion is 12.6 &%). 

Calculations have been performed within the same 
scheme as for FeBr 2 [3 and for NiO 0, OJ, [ll and 
other our earlier calculations for 4f/5f compounds [2J|. In 
Hamiltonian we take into account the crystal-field inter- 
actions (the octahedral crystal field with a small tetrag- 
onal distortion B 2 ), the spin-orbit interactions A s _ and 
inter-site spin-dependent interactions (ndd)- The result- 
ing low-energy electronic structure is shown in Fig. 1. 
In the magnetically-ordered state there appears a spin- 
gap of 466 K (=40 meV). The ordered magnetic moment 
amounts at T= K to 4.04 [is - it is built from the spin 
moment of 2.62 fis and the orbital moment of 1.42 \xb 
[2f| . The calculated orbital moment gives 35% contribu- 
tion to the total moment. The experimentally derived 
magnetic moment amounts to 4.0 fiB (26|. 

The tetragonal distortion is small and causes a slight 
splitting of the excited quartets. But this small distor- 
tion determines in our calculations the direction of the 
Co magnetic moment, in the CoO case with the c-axis 
compression, along the tetragonal direction. 

The formation of the magnetic state we have de- 
scribed for numerous compounds - let mention exemplary 
4f/3d/5f compounds ErNi 5 , FeBr 2 , NiO, FeO, UPd 2 Al 3 



and UGa2, results of which have been published start- 
ing from 1992. In all these cases the magnetic energy is 
much smaller than the overall CEF splitting. They are 
both ionic (FeBr 2 18], NiO, FeO [Tij) and intermetal- 
lic (ErNi 5 , UPd 2 Al 3 , UGa 2 ) [H compounds. A derived 
value of n= -43.8 T/fi B (= -29.4 K) in CoO yields the 
molecular field acting on the Co-ion moment at T = K 
as 176 T - it corresponds to the exchange field acting on 
the Co spin moment of 271 T and to the exchange energy 
of41.2meV. 

The calculated electronic structure very well describes 
the temperature dependence of the specific heat as is 
shown in Fig. 2, where the calculated data are com- 
pared with experimental data. The calculated temper- 
ature dependence of the specific heat of CoO contains 
three contributions: the 3d contribution c d (T) contain- 
ing both magnetic and Schottky contributions, the lattice 
and dilation contributions. The magnetic/Schottky con- 
tribution is calculated from the derived electronic struc- 
ture shown in Fig. 1 similarly to that discussed in Ref. 
[l8| for FeBr 2 . The dilation contribution is linear with 



temperature 27J - from thermal expansion and the bulk 
modulus we calculated its temperature coefficient as 3.4 
mJ/K 2 mol (we denote it as 7^ to be distinguished from 
the electronic Sommerfeld coefficient). The lattice con- 
tribution we approximate by the Debye function. We 
did not have intention to fit perfectly the data but using 
a value of the Debye temperature 9d of 600 K the ex- 
perimental data are nicely reproduced in the very wide 
temperature range, from to 1000 K. Here we give some 
values. At T=200 K c m =13.27, c tat =33.09 and c da 
=0.68 J/Kmol gives total heat of 47.02 J/Kmol, which 
should be compared with 46.03 HzJ (46.63 [H) J/K mol. 
At T=400 K c m =c Sc h =7.15, c lat =44.78 and c da =1.36 
J/Kmol gives total heat of 53.29 J/Kmol, compared to 
54.80 J/Kmol [H]. At T=600 K c m (=c Sch ) =5.35, Q at 
=47.62 and c da =2.04 J/Kmol gives total heat of 55.01 
J/Kmol, compared to 54.80 J/Kmol [13]. All of these val- 
ues are very close to experimental data of Refs [20|, [27| 
as seen in Fig. 2. 

The derived Debye temperature 9d of 600 K is quite 
similar to a value for 6b of 650 K we have obtained for 
NiO [l5j] - it is a very plausible result owing to the struc- 
tural and compositional similarity of both monoxides. 

We would like to note that all of the used by us param- 
eters (dominant octahedral CEF parameter lODq of 1.19 
eV (or i?4= -23 K), the spin-orbit coupling A s _ = -260 K, 
lattice distortions B 2 = -1 K) have clear physical mean- 
ing. The most important assumption is the existence of 
very strong correlations among 3d electrons what assure 
the preservation of the atomistic ionic integrity of the 
Co 2+ ion in the solid CoO. We stress the good reproduc- 
tion of experimental results like temperature dependence 
of the heat capacity with a A-type peak, the value of the 
magnetic moment and, in particular, its direction. 

The crystal-field interactions are relatively strong in 
CoO but not so strong to destroy the ionic integrity of 
the 3d electrons - it forms conditions for adequancy of the 
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QUASST theory. The calculated value lODq of 1.19 eV 
places CoO at value of 1.0 at the Tanabe-Sugano diagram 
drawn with Dq/B as abscissa [281. |2^|. 

III. CONCLUSIONS 

We have for the first time calculated low-energy elec- 
tronic structure both in paramagnetic and magnetic state 



as well as zero-temperature properties and thermody- 
namics. We consistently described magnetic properties 
of CoO in agreement with its insulating ground state. 
The orbital moment of 1.42 hb gives 35% contribution 
to the total moment of 4.04 hb at T =0 K. We have 
calculated from this low-energy electronic structure the 
temperature dependence of the specific heat being in nice 
agreement with experimental data. In our approach CoO 
is an insulator independently on distortions and the mag- 
netic order. 
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